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Abstract 

In this work, we design a nearly linear time discrete Morse theory 
based algorithm for computing homology groups of 2-manifolds, thereby 
establishing the fact that computing homology groups of 2-manifolds 
is remarkably easy. Unlike previous algorithms of similar flavor, our 
method works with coefficients from arbitrary abelian groups. Another 
advantage of our method lies in the fact that our algorithm actually 
elucidates the topological reason that makes computation on 2-manifolds 
easy. This is made possible owing to a new simple homotopy based 
construct that is referred to as expansion frames. To being with we 
obtain an optimal discrete gradient vector field using expansion frames. 

This is followed by a pseudo-linear time dynamic programming based 
computation of discrete Morse boundary operator. The efficient design 
of optimal gradient vector field followed by fast computation of boundary 
operator affords us near linearity in computation of homology groups. 
Moreover, we define a new criterion for nearly optimal Morse functions 
called pseudo-optimality. A Morse function is pseudo-optimal if we 
can obtain an optimal Morse function from it, simply by means of 
critical cell cancellations. Using expansion frames, we establish the 
surprising fact that an arbitrary discrete Morse function on 2-manifolds 
is pseudo-optimal. 

Classical Morse Theory cnun] analyzes the topology of the Riemannian 
manifolds by studying critical points of smooth functions defined on it. In 
the 90’s Robin Forman formulated a completely combinatorial analogue of 
Morse theory, now known as discrete Morse theory. The fact that Forman’s 
theory can be formulated in language of graph theory makes it possible 
to use powerful machinery from modern algorithmics to provide efficient 
algorithms with rigorous guarantees. It is worth noting that the reader can 
understand this work without any prior knowledge of Morse theory as long as 
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he understands the equivalent graph theory problem. Knowledge of discrete 
Morse theory is however useful for the more inclined reader who wishes to 
understand the context and wider range of applicability of this work. In 


subsection 1.2 we provide a quick overview of the graph theory setting of 


discrete Morse theory in order to enable the reader to make a quick foray 
into the core computer science problem at hand. 


1 Background and Preliminaries 

1.1 Discrete Morse theory 

Forman provides an extremely readable introduction to discrete Morse theory 

in m. 

Notation 1. The relation ’ is used to denote the following: t ~< a ^ t G 
a (fe dim r = dim a — 1. 

Notation 2 (The d-(d-l) level of Hasse graph). By the term, d-(d-l) level 
of Hasse graph % we mean the subset of edges of the Hasse graph that join 
d-dimensional cofaces to (d-1)-dimensional faces of Hasse graph. 

Definition 3. Boundary & Couboundary of a simplex a: We define 
the boundary and respectively coboundary of a simplex as bda = {r | r -< cr} 
cbd a = {p\(T < p} 

Definition 4. Discrete Morse Function: Let K, denote a finite regular 
cell complex and let C denote the set of cells of JC. A function T": T —)• M zs 
called a discrete Morse function (DMF) if it usually assigns higher values 
to higher dimensional cells, with at most one exception locally at each cell. 
Equivalently, a function F : C is a discrete Morse function if for every 
£ C we have: 

(A.) Afi{a) = ff{p G chda\F{p) < F{a)] < 1 
(B.) A/’ 2 ((t) = #{r £ bda \F{t) > T'(cr)} < 1 

A cell a is critical if J\fi{a) = A/ 2 (c) = 0; A non-critical cell is a regular 
cell. 

Definition 5 (Combinatorial Vector Field). A combinatorial vector field 
(DVF) V on C is a collection of pairs of cells {(a,/3)} such that {a"* -< 
yj(m+i)j each cell occurs in at most one such pair ofV. 

Definition 6 (Discrete Gradient Vector Field). A pair of cells {a^ -< 
F{a) > F{I3) determines a gradient pair. A discrete gradi¬ 
ent vector field (DGVF) V corresponding to a DMF F is a collection of 
cell pairs -< such that -< £ V ijf F(ft) < F{a). 
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Definition 7. We define V-path to be a cell sequence <7q , T( 


(m) (m+1) (m) 


^(m+1) 


cr, 


(m) (m+1) (m) 


0 


0 -; 


a. 


s.t. for i = 0,.. .q, {cjj -< rj G V, fXj r* 




CTj+i and Gi 7^ (Tj+i. T/ie V-path corresponding to a DMF F is a gradient 
path of F. 


Theorem 8 (Forman [6]). Let 1C be a CW Complex with a DMF T defined 
on it. Then JC is homotopy equivalent to a CW complex fl, such that Q 
has precisely one m-dimensional cell for every m-dimensional critical cell 
in 1C and no other cells besides these. Moreover, let Cm be the number of 
m-dimensional critical cells, fim the Betti Number w.r.t. some vector 
field V and n the maximum dimension of 1C. Then we have: 

The Weak Morse Inequalities: 


{A.) For every m € {0 .. .n}: we have Cm > fim- (1) 

(e.) CO - Cl ... + (-l)"c„ = /3o + (-1)^/3^ = x(/C) (2) 

The Strong Morse Inequalities: 

For every m G [0, n] : Cm - Cm-i • ■ ■ + (-l)”"co > fim - fim-i ■ • • + (-l)™'/3o 

(3) 

Notation 9. We shall denote the sum of Betti numbers by the symbol A and 
sum of number of critical cells by symbol T. In other words, 

n n 

2=0 2=0 


Notation 10. The symbol 0{n) is used to indicate nearly linear. It is given 
by 0{n) = n (log . 

Note 1.1. Civen a DCVF, we can use topological sort to obtain a total 
order on the cells and then assign (arbitrary) ascending function values 
to the sorted list of cells. This will give us a Morse function that agrees 
with the partial order imposed by the gradient vector field. Any such Morse 
function will have the same critical cells as the gradient vector field. Hence, 
we shall use the terms optimal Morse function and optimal gradient vector 
field interchangeably. 

Definition 11 (WMOC). Let T{N4) denote the sum of Morse numbers 
across all dimensions for the optimal DGVF on A4. We say that a fam¬ 
ily of simplicial complexes Q satisfies the weak Morse optimality condition 
(WMOC) when VAt G Q, T{M) = 0(1). In other words, T{M) |A1| 
uniformly \/AA G 12. 
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1.2 Graph Theoretic Reformulation 

Given a simplicial complex K., we construct its Hasse Graph representation 
'HfC (a-n undirected, multipartite graph) as follows: To every simplex a'^ £ JC 
associate a vertex G 'H/c- The dimension d of the simplex determines 
the vertex level of the vertex cr^ in T-Ltc. Every face incidence 
determines an undirected edge in T-L/c. Now orient the graph T-L/c 

to a form a new directed graph Ti/c- Initally all edges of Ti/c have default 
orientation. The default orientation is a directed edge —)• G Tiic 

that connects a k-dim node to a (k-l)-dim node Finally, associate 

a matching At to graph Ti/c- If an edge G M then, reverse 

the orientation of that edge to t cr'^ G Tiic- The matching induced 

reorientation needs to be such that the graph T-L/c is a Directed Acyclic Graph. 
A graph matching on T-Ljc that leaves the graph T-ijc acyclic in the manner 
prescribed above is known as Morse Matching. Table [T] provides a translating 
dictionary from simplicial complexes to their Hasse graphs. See [Figure 1 

1.3 Prior Work 

Joswig et al. |in| proved the NP-completess of the decision problem and 
posed the approximability of optimality of Morse gradient vector helds (for 
general dimensional complexes) as an open problem, by pointing out an error 
in Lewiner’s claim about inapproximability in m- Recently jl8j provided an 
0(log^ n) factor 0{n) time approximation algorithm for the optimal discrete 
gradient vector field (that minimizes the number of critical cells). Recently, 
Burton et al. [5] developed an FPT algorithm for optimizing Morse functions. 
Some of the notable works that seek optimality of Morse matchings by 
applying heuristics in general are m El 0 Eiiini ca m. The works that 
constitute more relevant prior work for us are those that achieve optimality 
by restricting the problem to 2-manifolds in nearly linear time muni and 
quadratic time respectively. Ours is however the first algorithm to compute 
homology groups of 2-manifolds with arbitrary coefficients in nearly linear 
time. 

2 Boundary Operator Computation 

The analytic formula for boundary operator is given in Forman|6|- The 
obvious interpretation of the formula gives an exponential time algorithm. 
We give an efficient 0{Kn) time algorithm where k is the total number of 
critical cells which is nearly-linear if the number of topologically interesting 
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DGVF on AT ~ Acyclic Matching on H 


{B, AB) (C, AC) (D, AD) {E, AE) 

{CD, ACD) {CE, ACE) {BE, ABE) {BC, ABC) 


Discrete Gradient Vector on complex K. 



Matching chosen is such that reoriented Hasse graph H remains acyclic 


Figure 1: Matching induced orientation of Hasse Graph 
Table 1: Graph Theoretic dictionary for Morse Matching 



Morse theory on Cell Complex K 

Graph theory on Hasse Graph Hic 

1. 

gradient Pair e V 

Matched pair of vertices {a,l3) € Hfc 

2. 

Dimension d 

Multipartite Graph Level d 

3. 

g.d-1 Y 

Default down-edge r ^ a 

4. 

g.d-1 5 ^.^ ^^d-1^ .j-d^g Y 

Matching up-edge a ^ r 

5. 

V-Path 

Directed Path 

6. 

Non-trivial Closed V-Path 

Directed Cycle 

7. 

CVF 

Matching on the Hasse Graph 

8. 

DGVF 

Morse Matching (i.e. Acyclic Matching) 

9. 

Critical Cell 

Unmatched Vertex C 

10. 

Regular Cell 

Matched Vertex ^ 
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features are small relative to the number of simplices. Hence, we have a 
pseudolinear time complexity algorithm for boundary operator Computation. 
We note of the following Theorem from Forman | 6 j: 

Theorem 12 (Boundary Operator Computation. Forman j 6 ]). Consider an 
oriented simplicial complex. Then for any critical (p+1 )-simplex fd set: 


Aid = ^ Paf} a 

critical ot(p) 

Pap= Y. 

■y&r{0,a) 

where r{ld, a) is the set of discrete gradient paths which go from a face in 
the boundary of jd to a The multiplicity 0{'y) of any gradient path 7 is equal 
toil depending on whether given 7 the orientation on fd induces the chosen 
orientation on a or the opposite orientation. The formula for the boundary 
operator above computes the homology of complex K. 

We observe that we need ’formal sums’ of critical cells at each critical cell. 
However, there is an advantage in calculating these formal sums for interme¬ 
diate regular cells as well since this can potentially speedup calculations at 
critical cells. Since topological sort also does ordering for us, we can start at 
the lowest valued critical cell. We proceed to the next higher valued cell and 
observe that we have two cases. 

Also, we assume that our complex has a pre-assigned orientation. The 
angular brackets <, > in the formulae above denote the pre-assigned orien¬ 
tation. Once the boundary operator is ready we use Smith normal form 
algorithm over a collapsed complex that is provably significantly smaller than 
the original complex, in a mathematically precise sense. 

Let us denote by Aa the boundary operator computation for cell u™'. We 
now make an inductive hypothesis that the computation of the A operator 
has been done for all the maximal faces / single coface (since they are all 
lower valued Morse cells). Then the value of the A operator for the new cell 
is calculated as follows: 

Case 1: All flow emanating from a cell goes out through its boundary 
faces. No lower-valued co-faces. 
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The first formula takes care of Case 1 where flow goes out through the 
faces of the boundary. Note that in the formula above, is a placeholder 

for non-critical faces (if any) of cr™', i.e. -< u™'}, which are not a part of 

the Discrete Gradient vector field which is equivalent to saying /(r) < /(cr). 
Similarly, is a representative for the critical faces (if any) ofu”*. This 

formula holds irrespective of whether cr™' itself is critical or non-critical. In 
case of computation of boundary of a critical fj™ such that m = 0 , i.e. when 
cr™ is a critical point, the boundary is null. 

Case 2: The cell has 1 lower valued coface. 


Act =< (9/3, cT > xA/3 (5) 

The second formula takes care of Case 2 when fj™ has a lower valued 
co-face /3™'+^. 

Case 3: The 0-dimensional cell a is the unique minima. 


Aa = 0 (6) 

Theorem 13 (Boundary Operator Computation: Correctness Proof). The 
Algorithm correctly computes boundary operator A. 

Proof. Note that, to begin with we start with a list of cells in an ascending 
total order. Let us call this list £. This total order is one of the total orders 
that is compatible with the partial order prescribed by the gradient vector 
field Y. If we assign the function value ’i’ i.e. the index of some cell C[i] 
to each cell in £, we essentially obtain a Morse function compatible with 
the gradient vector field. The first cell we process is one with the lowest 
function value (i.e. the unique minima). This cell is then followed by cells 
with increasingly higher Morse function values. To prove that the formulaic 
computation of the A operator as expressed in subroutine calcBdryOp() is, 
in fact, the same as expressed in Theoremwe proceed by induction. Let 
(Ti denote the unique minima. The base case of induction for Aui is trivial. 
Now suppose that for all cells in the set {ai,o' 2 , ■ ■ ■ u/}, we have correctly 
computed the boundary operator as prescribed in Theorem Now suppose 
we encounter cell ct/+i. Suppose that (T/+i has a lower valued coface ft i.e. 
((T /+1 -< /3 & ((T/+i,/3) G Y). Since ft has lower function value as compared 
to (T /+1 (by hypothesis), we conclude that ft = crj+i for some J < I. All 
paths emanating from (T/+i must go through ft. The orientation induced by 
some path 'yi ft ^ p from ft to some critical cell say p is i where t = ± 1 , then 

the orientation of path (T 7+1 p will be (5/3, (Tj+i) x l. Therefore, the total 
count of paths (with induced orientation accounted for) will be (5/3, a) x A/3. 
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Hence, the boundary operator computation done in calcBdryOp() is valid 
for the case when <t/+i has a lower valued coface. 

Finally, assume that cr/+i does not have any lower valued coface. Therefore, 
the flow leaving from (T/+i will be through each of its faces (except possibly 
one higher valued face). If it indeed has a (matched) higher valued face then 
flow will be entering it through that face and hence the face in question isn’t 
relevant in calculating the weighted sum of gradient paths that leave ct/+i. 
When consider lower valued faces of cr/+i, we make a distinction between 
faces that are non-critical and those those that are critical. If a face say 
cxj is critical, then clearly we are justified in directly including the entry 
OLj X {da, cxj) as part of our «formal sum» that makes up the cell boundary. 
As for the non-critical entries of the formula, namely [Atj x < da, Ti >], we 
impose an additional constraint (r*,^) G 'f'rn (as opposed to {£,,Ti) G T(„-i) 
in the summation. In doing so, we are ruling out all entries that would 
valid directed paths going out of ct/+i but those that won’t add up to make 
gradient paths as prescribed by Theorem [T^ Now since Ti is lower valued its 
boundary At* has already been calculated correctly by Induction Hypothesis. 
But clearly every gradient path emerging from uj+i must first pass through 
one of these r^’s. Also, for each of these gradient paths, the orientations will 
change precisely by the multiple of (tj,^). Therefore the weighted sum of 
(non-trivial) gradient paths from (T/+i will be the sum of all the contributions 
by boundaries of each of the non-critical faces r*. To complete the argument 
for the induction step, we note that these sums along with contributions 
from the critical faces of cr/+i takes into account each gradient path precisely 
once. Also, it is easy to see that multiplication by co-orientation at each step 
provides the weights to ensure that the final entry will decide the induced 
orientation. Hence proved. □ 


Theorem 14 (Complexity of Computing Boundary Operator). The com¬ 
plexity of computing the boundary operator is 0{T x M). 


Proof. For the Hasse graph H(V ,£) of a simplicial complex, £ <V xP where 
P is the maximum dimension of cells in the complex (which in our case is 
2). Therefore, |ifi| = 0(|V|). (It is easy to show that for a cubical complexes 
as well, number of edges is 0(|V|)).The complexity of computing topological 
sort of the oriented Hasse graph is 0(|V| -1- |<S|) which is same as 0(|V|), 
assuming that our input manifold is either simplicial or cubical. 

The for loop in Lines [8 21 of procedure calcBdryOp() costs at the most 
0{A) per iteration while the total number of iterations is 0(|V|). But 

A = T, the total cost of the for loop 


since m our case, using 


Theorem 29 


is 0(|V| X T). Therefore, complexity of computing boundary operator is 
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0(1 V| X T) = 0{N X T), since the number of vertices in the Hasse graph is 
same as number of cells in the complex (i.e. the size of the complex namely 

M). □ 


It is worth noting that in vast majority of the practical scenarios N ^ T, 
enough for us to assume that compared to the size of the complex, the 
’topological complexity’, T is nearly a constant. We therefore use the notation 
O(-) (where 0{T xN) = 0{N)) to indicate the nearly linear time complexity 
of boundary operator computation. 

3 Frames of Expansion 

3.1 Basic Formulation 

Notation 15. Let the set B{a) denote the 0-dim cells (vertices) and 1-dim. 
cells (edges) incident on a if a is 2-dimensional and let B(a) denote the 
0-dim cells incident on a if a is 1-dimensional. 

Definition 16 (Semigraph). A semigraph is a set of vertices and edges s.t. 
every edge may have either one or two vertices incident on it. 

Semigraphs generalize graphs in the sense that, in a graph, every edge is 
incident on precisely two vertices. 

Definition 17 (Frame of expansion of a critical cell). Given a critical cell 
a"' where (n>2), consider the set of all cells that can be reached from a, by 
following one of the gradient paths within the gradient vector field. We call 
this set the expansion set of critical cell a and denote it by a. The frame 
of expansion of a is the n — l-dim. boundary of a along with the n — 2 dim. 
cells incident on these boundary cells. We denote the frame of expansion of 
a by a. 

Definition 18 (Frame of expansion of a boundary cell). Given a regular 
boundary cell where (n>2), suppose that {f,,x) forms a gradient pair. 
Now consider the set of all cells that can be reached from by following one 
of the gradient paths within the gradient vector field. We call this set the 
expansion set of boundary cell f and denote it as The frame of expansion 
of f is the n — 1-dimensional boundary of We denote it as 

Note 3.1 (Method of addition of cells upon expansion). It must be 
noted that if there is an expansion along into cell , then we delete 
from the frame and the set B{w\t) is added into the frame. 
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Note 3.2. Suppose we are given a regular 2-cell vd^, s.t. the 1-cell 
G w‘^. The boundary of t namely B{t) consists of two vertices say A® 
and p^. Note that within the set B{w), there exist two non-intersecting 
paths that connect A and p. One path involves the singular edge t, the 
other path consists of edges belonging to the set B{w \ r) 


Definition 19 (connectedness, connecting path). Consider two cells i) 

in a complex /C”. We say that a and r are said to be Type 1 connected 
in complex JC if there exists a cell sequence 


(m) 

! Iq ! 


s.t. for i = 0 , 


1 . 7i 

■ q, (t>i < li h+i, 4>o = 


= r and (^i 7 ^ This sequence of cells, (pQ 


(m—1) 


is known 

as a connecting path. Analogously, we say that and are Type 2 
connected in complex KN if there exists a cell sequence 'y^\ (j)(^ 'y^\ 

<Pi , ■■■Tq , (Pq , 7g+i S.t. for I = 0,.. .q, Ji y (pi ^ 7 i+i, 70 = a, 

7q_l_i = r and 7* 7^ 7i+i- The sequence of cells, .. •7g™| is known as a 
connecting path. Finally, we say that, cr™- and are Type 3 con¬ 
nected if there exists a cell -< with Type 2 connectedness between 

rj and t. 

Finally, we say that a set of m — 1 and m dim. cells are said to form a 
connected set if for any pair of m — 1 dim. cells (alternatively, for any pair 
of m dim. cells) we can find sequence of connecting cells as prescribed above. 


3.2 Pseudocode for 0(n)-Time Algorithm for Computing Ho¬ 
mology of 2 -manifolds 

Notation 20. Given manifold A4, we use the notation to denote the 
d-dimensional cells of manifold N4. 

Definition 21 (Boundary faces, Coboundary faces). Given a complex 1C, 
if there exists a cell of dimension d s.t. there exists a unique (d + 1)- 
dimensional cell satisfying t) ^ w, then we call d a d-dimensional 

boundary face of complex 1C. Also, in this case, w is known as a {d -\- 1)- 
dimensional coboundary face of complex 1C. 

Definition 22 (Boundary and Coboundary). Given a complex 1C, the list of 
all d-dimensional boundary faces of 1C is known as the d-dimensional boundary 
of 1C. Also, list of all d-dimensional coboundary faces of 1C is known as the 
d-dimensional coboundary of 1C. 
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Definition 23 (n-flow). The set of gradient paths in vector field T' on 
manifold that involve alternating n-dim. and (n — l)-dim. cells is known 
as the n-flow of 'V 


Algorithm 1 Homology() 

1 : procedure CALcHoMOLOGY(At, A;) 

2: We use the mainFrame() subroutine to design a vector field Y on M.. 

3: We then use subroutine calcBdryOp() to calculate the boundary operator 
Ac for DGVF r. 

4: Finally, using chain complex implied by boundary operator Ac, we calculate 
homology of At (with coefficients coming from arbitrary abelian group A) 
using Smith Normal Form. 

5: end procedure 
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procedure CALcBDRYOp(At, A) 
topologicalSort('H, A, £, 'ASCENDING'); 

fji = A[l]; AcJi = 0; 
for 2 < f < |£|; cj := C[i] do 
if (fT,/3) is a gradient pair then 
Act =< 9/3, cr > xA/3; 

else 

Let Tj -< cr be the set of regular cells incident on a s.t. {Ti,a) ^ A; 
Let Oj -< O' be the set of critical cells incident on rr; 

Act = ^ Atj X < da, r* > + ^ cnj x < da, ai >; 

end if 

if cr is a critical cell then 
AcCJ := Act; 

end if 
end for 

Ac is the Morse boundary operator corresponding to vector field A. 

end procedure 


3.3 Frame Expansions: Correctness &; Complexity Proof 

Lemma 24. Suppose there exist two vertices and 7 ^ that are connected 
through edges that belong to some frame after a certain number of elemen¬ 
tary expansions. Then the two vertices will remain connected through edges 
belonging to that frame upon further expansions 
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Algorithm 2 Frame Flow 
1: procedure findCoBdry(A1'^) 

2; Let be the list of d-dimensional cells of manifold M.. Scan through the list 
If a cell has a face such that § is the sole coface of then 

add to = {ro G \ w has at least one boundary face. } 

3: return 
4: end procedure 

5: procedure ADDPAiRToVECTORFiELD(r, tD, 'y,M,B‘^, d) 

6: If r 7 ^ NIL and if r isn't already matched then do the following: 

7: (a.) Match r to-i?. (b.) Delete t? and r from and Al^ respectively. 

8: (c.) If r G .6^ then delete r from list B’^. (d.) Enqueue r in Q. 

9: (e.) Add (r?, r) to vector field A. 

10: end procedure 

11: procedure FRAMEFLOw(At^, d, A) 

12: Dequeue a cell w from If the dequeue operation with list B^ returns NIL 
then dequeue a cell from list At'^. 

13: repeat 
14: repeat 

15: if ro is a cell that has been dequeued from list B^ then 

16: if R is a boundary face of w then 

17: I Invoke addPairToVectorField() in order to add {v,w) to A. 

18; else 

19; I Delete vo from At'^. 

20: end if 

21: end if 

22; for each face di of w do 

23: If there exists t?* s.t. //j isn't part of any gradient pair of A 

24: Then invoke addPairToVectorField() to add {'di, fn) to A. 

25: end for 

26: Dequeue a cell from queue Q. Call it zu. 

27: until w 7 ^ NIL 

28: Dequeue a cell from queue B'^. Call it w. 

29: until w 7 ^ NIL 
30: end procedure 

31: procedure MAiNFRAME(At) 

32: Invoke findCoBdry() to find coboundary B^ of At^. 

33: Use frameFlow() to design vector field A on cells of At^. 

34; B[1 : numEars] := earDecompose(.y#) 

35; for 1 < f < numEars do 12 

36: Invoke findCoBdry(Bi) to find B^[i] 

37: Use frameFlow() to design vector field A on cells of ear Si. 

38: end for 
39: return A; 

40: end procedure 






Proof. By hypothesis, we assume that two vertices, say and 7 *^ are con¬ 
nected through edges belonging to the frame after a certain number of 
expansions. Therefore there exists a connecting path V connecting the two 
vertices. Suppose w.l.o.g., we expand along some edge into cell We 
have two cases. 

Case l:r^ ^ V. In this case, all edges in path V continue to belong to the 
frame after the expansion corresponding to gradient pair {t,w). Therefore, 
even after this expansion, a and 7 remain connected. 

Case 2:r^ G V. Suppose and are the vertices of r. Then there exists a 
path "Pi s.t. Pi C P connecting a and A. Also there exists another path V 2 
s.t. P 2 C P connecting p and 7 . However, from 3.2 we know that, A^ and p^ 
are connected through edges that belong to set B{w\t). The original path 
P consists of edges Pi U r U P 2 . Upon expansion, we have a new path namely 
Pi U {w \ r} U P 2 . Therefore, frame expansions maintain connectivity. □ 


Note 3.3 (2-Manifolds and Semi-graphs). A 2-manifold without bound¬ 
ary has the structure of a simple graph (unrelated to Hasse graphs) in 
the following sense: Let every 2-cell denote a vertex and let every 1-cell 
denote an edge connecting 2-cells. The manifold structure allows at most 
two incident 2-cells for every 1-cell, whereas not having a boundary im¬ 
plies that the incidence number is exactly two for every 1-cell. Now if we 
have a 2-manifold with boundary, then the boundary 1-cells will have only 
one incident 2-cell whereas all other 1-cells will have two incident 2-cells. 
Therefore a 2-manifold with boundary has the structure of a semi-graph. 
For a given 2-manifold A4, let us denote the semigraph structure by Qs{Ai). 


Lemma 25. Every vertex belonging to manifold Ai is included in the frame, 
when all expansions are processed. 


Proof. From Note 3.3 we know that the 2-cells and 1-cells of a given 
2 -manifold Ai forms a semi-graph structure which we denote by Qs{Ai). 
We use the following convention: If a 2-cell, say r is included in some gra¬ 
dient pair belonging to vector held P or if r is the start cell of procedure 
frameFlow() described in Algorithm 2 the we say that vertex r is traversed. 


Case 1: Suppose Ai has no 1-dim. boundary faces. Then Qg[Ai) as¬ 
sumes the structure of a connected simple graph. In this case, the procedure 
frameFlow() described in Algorithm 2 we begin with some starting 2-cell w. 
While scanning through all the faces of w, if we hnd a face D s.t. p ^ w and 
p >- {1 and p isn’t part of any gradient pair, then we traverse p by adding 
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gradient pair i?, // to vector field 'V and add /r to a queue. Having processed 
all faces of tu, we dequeue a cell, say ronew from the the queue. We process 
in exactly the same way as we process w. And we keep doing this till the 


UJ I 


queue is empty. Clearly, this is equivalent to a breadth first traversal on graph 
Gsi-M.). Given the fact that all vertices of a graph are traversed in a breadth 
first traversal, we conclude that except for the start cell, all other 2-cells are 
part of some gradient pair. When the start cell w is added, the expansion 
frame consists of B{w). Every time we add a gradient pair ("i?,/u) to , we 
delete from the frame and add B{fi \ i9) to the frame. Since every vertex Vi 
belonging to M is part of H(/r \ •d) for some 2-cell /r, we see that each vertex 
Vi becomes part of the expansion frame at some stage of the construction of 
the frame. When new gradient pairs are processed, we may delete 1-cells from 
out frame, but 0-cells are never deleted. So, all vertices of A4 eventually be¬ 
come part of the expansion frame. See Figure 11 and Figure l^for an example. 


Case 2: Suppose Ai has some 1-dim. boundary faces. In this case, Qg[Ai) 
has the structure of a possibly disconnected semigraph. If the manifold 
has a coboundary face, say Bi then for the first connected component of 


Qs{Ai) in lines 16 17 of Procedure frameflowf) in Algorithm 2 we add the 
boundary-coboundary pair to vector field A. Following that, the cells that 
are in the same connected component are added to the vector field in a 
manner similar to Case 1. If there exists another connected component, then 
surely such a connected component must have at least one coboundary face. 


In lines 29 and 28 of Procedure frameflowf) in Algorithm 2, we check if such 


a coboundary face exists. If it does exist then in the loop 13 29 we process 
the every connected component in the same way as we process the very first 
one. Given the fact that all connected components of semigraph QflAi) are 
processed, every 2-cell Wnew in each of these components is part of vector 
field Y. Suppose Wnew is paired with some 1-cells 7?new each time then, each 
time we delete 7?new from the the frame and add B{p\ 'dnew) to the frame. 
When new gradient pairs are processed, we may delete 1-cells from out frame, 
but 0-cells are never deleted. Since every vertex is incident on at least one 
of the 2-cells in one of the connected components, we establish the fact that 


all vertices eventually become part of the expansion frame. See Figure 8 
Figure 9 and Figure 9] for such an example. 


□ 


^Example: For a connected manifold K., this may happen, for instance when, say one 
the 2-cells A is connected to other 2-cells only by the medium of 0-cells while the 1-cells of 
A are not shared with other cells. See [Figure 8| [Figure 9| and [Figure 9| for another such 
example. 
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Note 3.4. Every connected component ofQs{-M) has a 2-cell which shares 
a 0-cell with a 2-cell from another connected component. If we imagine 
every connected component ofQs{M) as a vertex and every shared 0-cell as 
a hyperedge then we get a connected hypergraph that we write as Hd-M.)- 
We know that HciM) is connected because if this were not the case then 
clearly At itself will have more than one connected components. We call 
HciM) the component hypergraph of M. 

Consider a gradient vector field V assigned to a manifold At. First consider 
the case when At is a manifold without boundary. Consider a critical cell cd. 
Note that before we do any expansions, B{ol) is our original S. Let a t < ol. 
If (r, O') G V then we can consider it as an expansion along r to the cell a. 
Now, as per the definition of frame expansion, we add the set {B{a) \ {r}} 
to S and we delete {r} from a. Therefore, 

3 = a-{r} + {fi(o-)\{r}} 

But, this is same as saying, a = {B(a) \ {r}} + {B{a) \ {r}}. Clearly, the 
sets {B{a) \ {r}} and {B(a) \ {r}} are themselves connected and both these 
sets have a common boundary namely B{t) (the boundary of r). Therefore, 
expansion along r preserves the connectivity of S. Now, given our intermediate 
stage a, if any of the 1 dim. cells say OjGa forms a gradient pair with a 2 
dim. cell wf then by expansion we have, 

a = a - {di} + {B{wi) \ 

Each time we observe that the boundary of a and the boundary of {B{zui) \ 
{t?i}} is, in fact, the same as the boundary of Oi namely B{'di). Therefore, 
upon expanding the frame a along Di, connectivity of a is preserved and 
a continues to be a 1-manifold without boundary. Note that owing to the 
manifold nature of AI, {B{zui) \ {d*}} never contains a (l)-dimensional face, 
say Oj (where j < i), along which a was previously expanded. Therefore, 
because AI is a manifold, the two encounters of dj can happen in two different 
ways, namely: 

Case 1: While constructing a through expansions, any face "dj can be 
encountered at most twice - once when it is included in a as part of some 
{B{zuk)} {k < j) and a second time if and when we expand along Even 
as we expand along -dj, the two vertices of "dj stay connected. 

Case 2: The other possibility of two encounters for the face -dj is when it 
is included in a as part of some {B{wk)} {k < j) and some {B{'[Uh)} {h < j). 
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In this case, we never expand along 

If Ai has boundary then our start cell is a coboundary face and upon first 
expansion, the frame is a manifold with a boundary. Applying the reasoning 
above, for a given connected component of Qs{Ai), the frame of expansion 
restricted to a single connected component of Qs[Ai) is a connected 1-manifold 
without boundary. To arrive at the more general conclusion that the frames 
of expansion of all connected components of Qs{Ai), pieced together form 
a single connected 1-complex connecting all 0-cells of manifold AI, we have 
the lemma below: 


Lemma 26. Given any sequence of elementary expansions, the frame of a 
critical cell of a manifold Ai is always a connected set. Following the final 
expansion, the frame consists of a set of edges that connects all vertices of 
the complex. 


Proof. Consider without loss of generality, that Ad is a manifold without 
boundary. Then Qg(^A4) has a single connected component. Every ver¬ 
tex within the frame that was previously connected, stays connected by 
[Lemma'^ Since, for a manifold without boundary, the frame always has 
a single connnected component at every stage of expansion, and since by 
[Lemma ^ all vertices become part of the frame, we arrive at the conclusion 
that all vertices of the frame form a single connected component at the 
conclusion of all expansions. 

The other case, when Qs{-M) has several connected components, we first 
observe that the frames of each of the connected component stays connected 
by the same logic as in case of expansions of manifolds without boundary. 
Also, we observe that in such cases, every connected component will have 
a 2-cell which is connected to another connected component via a common 
0-cell. In fact, if there exist vertices Va and Vf, in two different connected 
components Ca and Cb- Ca and Cb may be interpreted as vertices in the 
hypergraph then we can first determine a path between Ca and Cb within 
the component hypergraph Hc(AI). Now, every vertex Ci in the path is a 
connected component and every hyperedge is a shared 0-ceII Vi. If the path 
is written as Cq, fi, Ci,.. ■ ,Vi, Ci,... Vn, Cn where Cq = Ca and Cn = Cb- 
Then for every C* 1 < z < n, we can determine an internal path (part 
of the expansion frame) in graph Qs{Ai) between Vi and Uj+i. Finally, in 
graph Qs{Ai), we can find a path between Va and vi within component Cq 
and a path between Vn and Vb within component Cn as parts of expansion 
frames within those components. If we piece together each of the paths from 
expansion frames of various components of Qs{Ai) along the path in the 
hypergraph %c{Ai), we get a path connecting any two vertices Va and Vb 
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such that every edge in the path is part of the expansion frame. From this 
we conclude that all vertices in the complex are connected to each other 
through edges that lie entirely in the expansion frame. In other words, the 
expansion frame is a single connected component that connects all vertices 
of the complex. □ 

Lemma 27. Applying the frame based algorithm on a 2-manifold gives us: 

C2 = /32 

Proof. Casel: ,02 = 1 Suppose the 2-manifold does not have a boundary. 
Then clearly /32 = 1. Now we will prove that in this case, C 2 also equals 1. 
Recall that Qs{Ai) takes the structure of a simple connected graph and the 
procedure framellow is equivalent to a breadth first traversal that begins with 
a start cell tu, where w is not included in any of the gradient pairs. However, 
subsequently every neighboring 2 -cell is paired with a 1 -cell and added to a 
queue. The neighbors of the dequeued cell are then scanned and if unpaired, 
they are paired with the connecting 1-cell as before. This process is continued 
till all 2 -cells are exhausted (which happens at the conclusion of the breadth 
first traversal). Hence all 2-cells except the start 2-cell varpi form a gradient 
pair with some 1 -cell, giving us C 2 = 1 . 

Case2: /02 = 0 Now, consider the case when the 2-manifold has a boundary. 
So, we have /02 = 0 and we will prove that C 2 also equals 0. Note that, in 
this case, Qg[AA) has one or more connected components s.t. each of the 
connected components has at least one coboundary face. For every component 
a coboundary face is selected as a start cell and paired with a boundary face 
to give a gradient pair. Subsequently, as before neighboring 2-cells are paired 
with connecting 1-cells if they haven’t been paired before. Newly paired 
2 -cells are queued and this process continues till all 2 -cells of the connected 
component are exhausted. In other words, every 2-cell of every connected 
component is part of a gradient pair giving us C 2 =0. Hence proved. □ 

Note 3.5. Let be the coboundary of residual complex Ai. B^[i\ is part 
of the coboundary that intersects with ear Si. i.e. B^[i\ = B^ H Si. 

Lemma 28. If the complex is made up of a single connected component, then 
the frame based algorithm gives its cq = 

^The case when the complex is made up of several connected components can easily be 
dealt with by applying the algorithm independently to each of the components. In that 
case Co = £ where £ is the number of connected components 


17 



Proof. From Lemma 26 we know that the frame of expansion consists of a 
single connected component that connects all 0-cells in the manifold. This 
frame is divided into ,yV several ears say Tj. Every ear is a 1-dimensional 
manifold. Suppose that we have an open ear then we have 1-dimensional 
coboundary face in such a ear which we pair with a 0-dimensional boundary 
face. Subsequently, we follow a path which matches the incident unpaired 
0 -cell to a neighboring 1-cell and we keep doing this until all 1-cells of the ear 
are exhausted. Now suppose that we have a closed ear. Then we remove one 
of the 1-cells from the ear (i.e. make it critical). This disconnects the ear 
into two connected components. We treat these two components of the ears 
as separate and proceed as in case of open ears. We now make an inductive 
argument to prove that the first ear leaves a critical 0-cell. Subsequent 
addition of ears do not add any criticalities. To see this consider the base 
case in which we design the flow for the first ear. Here, the flow stops when 
all 1-cells are exhausted. In this case, for the final 1-cell /i, there is one 0-cell 
which gets paired with /x and another incident 0-cells which remains unpaired. 
It is this 0-cell that becomes the sole critical 0-cell. For induction consider 
the inductive hypothesis that k-ears have been attached and the number of 
critical cells remains 1. Now suppose that the (k+l)th ear is attached. If 
the (k+l)th ear is open then the flow stops with a 1-cell on which one of 
the incident 0-cells Vi belongs to a ear Tj where i < {k + 1). Either Vi is 
the sole critical 0-cell or it is paired to another 1-cell belonging to Si (by 
inductive hypothesis). Now, suppose that the (k+l)th ear is closed. Then 
having detached a 1-cell (which is made critical), we have two disconnected 
components. For each of the connected components, the flow emanating from 
subsequent pairing of 0-cells to 1-cells stops when a 1-cell is incident on a 
0-cell Vj belonging to a ear Sj where j < (A: -|- 1). Once again by inductive 
hypothesis either Vj is the sole critical 0-cell or it is paired to another 1-cell 
belonging to Sj. From this we conclude that cq = 1 on attachment of all 
ears. □ 


Theorem 29. For the frame-based vector field design algorithm, each Morse 
number equals the Betti number, i.e. 


Ci = Pi 

Proof. From Lemma 27 and Lemma 28 

we have 

C 2 = P 2 and Co = Po 

respectively. Now, using Equation 2 

in 

Theorem 8 

we have ci = Pi. Thus 


we have Cj = jdi for all i. 


□ 
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3.4 Discussion on Complexity 

Finding coboundary of can be found in linear time by going through all 
2-cells in A4^. Finding coboundary of ears of ^ can be found in constant 
time by mainting a proper data structure. The ear decomposition of residual 
complex ^ (which has the structure of a graph) itself takes linear time. 

Adding a gradient pair to a vector held takes constant time. The queueing, 
dequeueing and deletion operations also can be done in constant time by 
maintaining appropriate data structures. 

The only nontrivial procedure in the algorithm is frameflow(). Now the 
frameflow() procedure can be construed as breadth hrst traversal on a semi¬ 
graph. We apply this procedure once on A4^ and once on each of the ears 
of When traversals from all ears are counted, we observe that every 
edge of ^ is encountered only once and every vertex Vi is encountered T){vi) 
number of times where T){ ) indicates degree of a vertex. So, if we sum over 
all vertices and edges, the total complexity of frameflow() when applied over 
^ is linear in the number of edges of ^. Hence, we see that the design 
of optimal discrete gradient vector held using expansion frames takes linear 
time. 

4 Pseudo-optimality of Random Morse functions 

In this section, we establish the surprising potency of critical cell cancellations 
in case of 2-manifolds by using frames. 

Definition 30 (Pseudo-optimal Vector Field). We define a DGVF to he 
pseudo-optimal if the optimal DGVF can he obtained from it merely via 
critical cell cancellations. 

Definition 31 (Stable, Unstable Manifolds). The stable manifold of a critical 
cell are all the non-critical cells of dimension q and q 1 with gradient 
paths ending at a'^. The unstable manifold of a critical cell are all the 
non-critical cells of dimension q and q — 1 with gradient paths starting at 
and ending at that particular non-critical cell. 
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Algorithm 3 Optimal DGVF Redesign Using Critical Cell Cancellations 
1 : procedure KiNGREv(j^9, At, C,y, O') 

2 : repeat 

3: Suppose critical cells and have gradient paths to/from saddle 7 ^. 

4: Subroutine sharedSaddle() finds such a pair { 7 ,cj} for given M. 

5: If 7 7 ^ NIL, then cancel critical pair ( 7 ,( 7 ) 

6: until (7 7 ^ NIL) 

7: If g = 2 AND has a unique path to then cancel critical pair 

8: end procedure 

9: procedure FixBDRY((i,p, At,, 6 ^) 

10 : for 1 < i < \I3‘^\ do 

11: Let Bi := and let bij be a boundary face of Bi. 

12 : if {bij,Bi) ^ V AND bij is critical then 

13: Let {9i, Bi) be a gradient pair 

14: if / NIL and 6i is not a boundary face of Bi then 

15: I Find a gradient path from some critical cell to 9i and reverse it. 

16: end if 

17: Add gradient pair {bij,Bi) to vector field A 

18: end if 

19: end for 

20 : end procedure 

21: procedure FiNDKiNG((i,p. At, 6'?, i) 

22 : if(p = 2 AND 7 ^ NIL) OR (p = 1 AND i = 1) then remFrom(C^, .^); 
23: else if p = 1 then ^ := 6'?[1]; remFrom( 6 '?, .^); 

24: else ^ ;= NIL; 

25: end if 

26: := selectRandomly(CP); return A; 

27: end procedure 

28: procedure PROGESsCoMPLEx(At, i, (i,p) 

29: C[1 : d] ^ identifyCritical(At, A) 

30: findBdryO finds b^~^ the cobdry. and bdry. of At resp. 

31: fixBdry((i,p, At,C‘^,,B'^) 

32: while (.^ = findKing(d,p. At, C^, b'^~^,i) 7 ^ NIL do 
33: kingRev(Jt, M,c,r,p) 

34: end while 
35: end procedure 

36: procedure KiNGFLOw(5[)t,d, A) 

37: Divide sot into manifolds Ati, At 2 , 20 -, At if s.t. At* n At^ is 0-dimensional. 
38: for 1 < i < AT do processComplex(Ati, i, 2, 2) end for 
39: iS[l : numEars] := earDecompose(.y#) 

40: for 1 < f < numEars do processComplex(<S[i], f, 1,0) end for 
41: end procedure 








Note 4.1. Given a connected pseudomanifold complex divide 
into several connected components A4 i,j\ 42, ■ ■ ■ ,-M.K s.t. Mi Ci Mj is 
0-dimensional, i.e. any of the two manifolds (with boundary) Mi,Mj 
may intersect only along points (but not along edges). IfHHisa manifold 
without boundary then will have only one connected component. Mi 
are essentially the connected components of the semigraph Qg[M) defined 

m liVote 3.3[ 


Lemma 32. If iJJl is a manifold without boundary then after invoking the 
procedure kingRev(), we obtain a connected expansion frame. Moreover, 
C2= h = 1- 

Proof. Suppose a vector field Y on manifold 501 (without boundary) has a 


single critical cell. Then from Lemma 26 we get a single connected expansion 
frame connecting all vertices of Instead, if 901 is a manifold without 
boundary and if we have more than one critical 1 -cells, then consider the 
unstable manifold of some chosen critical cell Since the unstable manifold 
of ^ doesnot include the entire manifold 911, the stable manifold has a 1- 
dimensional manifold as its boundary. From [ 12 ], we know that, if M is an 
n-dimensional manifold with boundary, then the boundary of M is an (n- 1 )- 
dimensional manifold (without boundary) when endowed with the subspace 
topology. Therefore, the boundary of the unstable manifold is a 1-dimensional 
manifold without boundary (i.e. it consists of one or more disjoint circles). 
Clearly the 1-cells belonging to this boundary are not part of the 2-flow of 
Y, else they wouldn’t be part of the boundary of the unstable manifold of K. 
So, the 1-cells belonging to this boundary are either part of the 1 -flow of Y 
or they are critical. Consider one of the disjoint circles that forms part of 
the boundary of the unstable manifold. If all the cells on this circle are part 
of the 1-flow then it will form a cycle. Hence there exists at least one critical 
1-cell on the boundary of the unstable manifold. Let 7 be a critical that lies 
on the boundary of the unstable manifold of .h. Clearly, there exists only one 
gradient path from to 7 . 7 is also incident on a 2 -cell say ui that does 
not he in the unstable manifold of R. Suppose ai is itself a critical 2-cell, 
then 7 lies on the boundary of unstable manifolds of the two critical 2 -cells 
^ and (Ti. Otherwise suppose that ai is matched. Because the simplicial 
complex 911 is a manifold, it is possible to trace any inverted gradient path 
on 911 (such a unique inverse gradient path exists). Therefore, we trace the 
inverted gradient path 7 , cri,... until we reach a critical 2 -cell (say ak) from 
which this path emanates. In any case, we can find a critical 1-cell 7 which 
is shared by critical cells ^ and some other critical 2-cell say cr. In this case. 
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because gradient path from u to 7 is unique we can invert this gradient path 
as shown in Line Definition 6 of Procedure kingRev() of Algorithm 3 Once 


this cancellation is done, the unstable manifold of a becomes part of the 
new unstable manifold of Once again we search a critical 1-cell 72 on the 
boundary of the unstable manifold s.t. which also lies on the boundary of 
unstable manifold of some other critical 2-cell (distinct from ^). If such a 
pair of critical cells is found then we cancel it and this procedure is repeated 
until all critical 2 -cells belong to the unstable manifold of ^ (or alternatively 
all critical 1-cells have two gradient paths from R.) Basically this means 
that is a manifold without boundary that has a unique critical 2 -cell. i.e. 
C 2 = 1. Since, HR is a 2-manifold without boundary, /32 = 1- Finally, from 
Case 1 of Lemma 27, we arrive at the conclusion that the expansion frame is 
a connected 1-manifold that includes all 0-cells of HR. □ 


If HR is a manifold without boundary then Gs{-M.) has a single connected 
component and the for loop described in Lines 38||M of Procedure kingFlow() 
in [Algorithm 3 gets executed only once. Also the while loop described in 
of Procedure processComplex() in Algorithm 3 gets executed 
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Lines 

only once for manifolds without boundary. This is because for any critical 
2 -cell .ft, you always find another critical 2 -cell a s.t. both ft and a have a 
gradient path to a common 1 -cell 7 unless the unstable manifold of ft covers 
the entire manifold HR. 

The situation is however much different for a manifold with boundary. For 
such a manifold the for loop and the while loop may run several iterations. 


Lemma 33. // HR is a manifold with boundary then after invoking the 
procedure kingRev(), we obtain a connected expansion frame. Moreover, 
C2 = P 2 = 0 . 


Proof. We will examine the effect of the algorithm on one of the connected 
components Mi of QsiM). Consider the unstable manifold of a critical 
2-cell ft. From | 12 | . we know that, if M is an n-dimensional manifold with 
boundary, then the boundary of M is an (n-l)-dimensional manifold (without 
boundary) when endowed with the subspace topology. Hence, the boundary 
of this unstable manifold will be a 1 -manifold without boundary (i.e. a 
disjoint set of circles). The 1-cells on any one of these circle are involved only 
in 1-flows or they are critical. But all, cells of a circle can not be involved in 
1 -flow as this would lead to a cycle in the vector held. So, every circle must 
contain a 1 -cell, say 7 that is critical, ft has only one gradient path to 7 . 
There exists a second gradient path that ends at 7 . This gradient path either 
emanates from another critical 2 -cell say a or it emanates from a boundary 
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face. Assume the case where a path to 7 emanates from a. In this case, the 
pair ( 7 , (t) is detected and cancelled in Lines [3]j^ of Procedure kingRev() 


m 


Algorithm 3 In fact, every such pair ( 7 , cr) for a given is detected and 

So 


cancelled in the loop Lines 3]|6 of Procedure kingRev() in Algorithm 3 


finally every critical 1 -cell say </> in the boundary of the unstable manifold 
of ^ will have a second path emanating from a boundary face. In this case, 
the pair of critical cells (</>, is detected and cancelled as shown in Line 
[^of Procedure kingRev() in Algorithm 3 Suppose that Aii continues to 


have critical cells that are not cancelled, then a new critical king cell M is 
selected and the same procedure as described above is repeated in a loop 
shown in Lines 34 of Procedure processComplex in Algorithm 3 We 
exit from the loop provided there are no other critical 2 -cells to process in the 
list CP. In case of manifolds with boundary every critical 2-cell processed as a 
king cell ^ is itself cancelled along with cancelling all the neighboring critical 
2-cells that share gradient paths to the same saddles as Having processed 
Mi in this manner, we are assured that eventually Mi has no critical 2 -cells. 
In fact every 2-flow for Mi emanates strictly from boundary faces. Using an 
argument similar to that in Case 2 of [Lemma 27| we know that the frame 
of expansion of a boundary face is a connected set. Consider the first such 
boundary face 61 , with a frame of expansion which is a connected 1 -manifold. 
Every 1-cell belonging to the frame of expansion of bi has a second gradient 
path emanating from other boundary cells {bi}. Since the Mi is a manifold 
with boundary, given any pair of boundary faces bi, bj, we can find a type 2 
connected 2-path between them. Consider all the 1-cells in some such type 
2-connected path between bi and bj. Every 1-cell either lies in the frame 
of expansion of two boundary faces or is involved in 2 -flow with a regular 
2-cell. This gives us a sequence of frames of expansion of boundary faces 
bki,bk 2 , ■ ■ ■, bk^ that are sequentially pairwise connected and s.t. bk^ = bi 
and bk^ = bj. Since this procedure can be applied to any two boundary faces 
(with expansion frames), we conclude that the set of frames of expansion of all 
boundary faces is a connected set, which we refer to as the expansion frame 
oi Mi. To see that the frame of expansions of all Mi form a single connected 
set, we consider the component hypergraph T-Lc{M). We then use the same 
line of reasoning as used in Lemma 26\ to conclude that the expansion frame 
of 971 is connected. Also, following all critical cell cancellations since there 
are no more critical 2 -cells for Mi, we have C 2 = ^2 = ^ for each Mi. So, we 
have also have C 2 = P 2 = 0 for 971 □ 
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Note 4.2. The residual 1-complex ^ is essentially the expansion frame 
of iXfl following cancellation of critical cell pairs of dimensions 1,2. Since 
there exists a preordained 1-flow (without cycles) on ^, clearly given the 
mechanism of discrete Morse theory, there must exist at least one critical 
0-cell. (A sub-optimal 1-flow may have more than one critical 0-cells. But 
at least one is guaranteed.) The first ear is chose to be one that includes 
at least one of these critical 0-cells. Also the ear decomposition follows a 
special procedure. The number of ears are determined by the number of 
unstable manifolds of boundary 0-cells and critical 1-cells. The first ear is 
either an unstable manifold of a boundary 0-cell or a critical 1-cell that 
includes at least one critical 0-cell. The second ear is a 1-manifold that is 
incident on a 0-cell that belongs to the 1^* ear and includes all the 0-cells 
and 1-cells of an unstable manifold of a boundary 0-cell or a critical 1-cell 
that aren’t already included in the first ear. The ear is a 1-manifold 
that is incident on a 0-cell that belongs to one of the previous {k — 1) ears 
and includes all 0-cells and 1-cells of an unstable manifold of a boundary 
0-cell or a critical 1-cell that aren’t already included as part of the previous 
(k — 1) ears. Every ear (apart from the first ear), has at least one 0-cell 
in its 0-dim. boundary whereas every ear may have at most two 0-cell in 
its 0-dim. boundary. The first boundary cell of the ear 6^[1] is incident on 
one of the previous ears. The second boundary cell [2] may or may not 
be incident on any of the previous ears. 


Lemma 34. On applying a series of critical cell cancellations, the connected 
expansion frame has cq = flo = 0 


Proof. We shall make an inductive argument. The idea is that the first ear 
will have a 0-cell that is critical. Subsequent ears attached to the first ear 
have no 0-dimensional critical cells. Note that all ears are 1-dimensional 
manifolds (topological circles or topological line segments) 

Base Case: Suppose that we start with the first ear. Suppose that the first 
ear is a closed loop (i.e. a topological circle). From Note 4.2 our first ear 
has at least one critical 0-cell. Suppose we call it In this case g becomes 
our king critical cell K. If there exist two gradient paths to g from a saddle, 
then clearly we do not have any criterion for cancellation. Instead if we have 
a single gradient path from the saddle 7 ^ and suppose there exists another 
gradient path from the 7 to some other minima then from Lines 3]|6 of 
Procedure kingRev() in Algorithm 3 we cancel critical pair ((^, 7 ) and as a 
result have a single critical 0-cell in the first ear. The last possibility the first 
ear consists of the unstable manifold of a critical 0 -cell ? that g^ 
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Induction step: By the inductive hypothesis, we have processed {k — 1) 
ears so far and for all the {k — 1 ) ears taken together, we have only one critical 
0-cell (namely the one that was encountered in the very first ear.) Now, we 
need to establish that on attachment of the k^^ ear we do not introduce any 
new critical 0-cells. Note that for k^^ ear we start with 6°[1] as the king 
cell .^, where is incident on one of the previous ears (i.e. it may either 
be our original critical cell or it may be some regular 0 -cell from one 
of the earlier (/c — 1) ears. Like all other ears, the k^^ ear is an unstable 
manifold of a boundary 0-cell or a critical 1-cell. If it is the unstable manifold 
of a boundary 0 -cell <;■ then we do not have anything to prove as the flow 
for this cell will simply start with and end at without introducing 
any criticalities. If the k^^ ear is topologically a loop, then has two 
gradient paths from some saddle 7 and hence the criterion for cancellation is 
not satished. Yet another case is when and 6 *^ [2] are both incident on 
one of the earlier k ears. In this case, is either or a regular 0-cell and 
6 ^ [ 2 ] is certainly a regular cell. Also, there does not exist any other critical 
0 -cell in this ear because the k^^ ear, in this case, is an unstable manifold 
of a saddle. The only interesting case is when k^^ ear is topologically a line 
segment s.t. [ 2 ] is critical and the saddle 7 has one gradient paths each to 
6 ^[ 1 ] and 6 ^ [2]. Since, for k^^ ear we start with 6 ®[ 1 ] as king cell .ft, we end up 
cancelling 7 and [ 2 ], making the k^^ ear an unstable manifold of boundary 
cell 6 ® [2]. In each of the cases, we ensure that either the k^^ ear did not have 
any critical 0 -cell to begin with or if there does exist a critical 0 -cell, then it 
is cancelled. Hence proved. □ 


Theorem 35. Every discrete gradient vector field on a 2-nianifold is pseudo- 
optimal. 


Proof. Suppose that at the end of the hrst call to Procedure findKingO 
from Line 
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of Procedure processComplex() in Algorithm 3 is not 
NIL. Then, we claim that the unstable manifold of .ft^ does not have any 
critical 1-cells that are boundary faces. This is because, if did have any 
boundary critical 1 -cells in its unstable manifold, it would have got cancelled 
in the loop shown in Lines lOpT in Procedure fixBdry() in Algorithm 3 In 
fact, more generally every critical 2-cell at the end of first call to Procedure 
findKingO will have no critical boundary 1 -cells in their respective unstable 
manifolds. If 501 is a manifold without boundary, then by [Lemma 32 


we 


have C 2 = /32 = 1 and we get a connected frame of expansion in form of 
residual complex ./#. Instead, if 501 is a manifold with boundary, then from 


Lemma 33 we obtain C 2 = /?2 = 0 and a connected frame of expansion in form 


of residual complex Given a connected frame of expansion guarantees 
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that we have cq = /3o = 1- Finally, using Weak Morse Inequlity we obtain 
Cl = /3i. Hence, we prove that for a 2-nianifold, given an arbitrary vector 
field 'Vx merely by using critical cell cancellations, we may obtain the optimal 
vector field In other words, every gradient vector field on a 2-manifold is 
pseudo-optimal. □ 


The Topological Explanation for simplicity of com¬ 
putation of i7(At^,A) 


We compute homology using [Algorithm 1 Also we assume Weak Morse 


Optimality Condition as defined in [Definition 11 on the input. 


As we can see from arguments in section 3 the topological explanation 
for simplicity of computation of homology groups for 2-manifolds is: 

1. On 2-manifolds optimal Morse functions are perfect. In fact, 2-manifolds 
admit readily computable perfect Morse functions. 

2. A 2-manifold has optimal C 2 = 0 or C 2 = 1, which can be figured out in 
linear time by examining whether or not it has a boundary. 

3. We define and apply frames of expansion an elementary homotopy 
theory construct to design our algorithm. 

It can be seen that irrespective of what traversal method we use to 
traverse the graph like connectivity structure of 1-cells and 2-cells of a 
2-manifold, the frame of expansion remains connected. Furthermore, 
this connectivity guarantees that optimal cq = 1. 

Finally weak Morse Inequality guarantees that our ci is optimal, i.e. 
Cl = /3i. 

Moreover, our dynamic programming based boundary operator compu¬ 
tation algorithm is pseudo-linear time (which becomes strictly linear 
assuming WMOC). 

7. Finally, assuming WMOC, the application of Smith Normal Form (a 
supercubical time algorithm) on input of constant size is inexpensive. 
The pseudo-optimality of arbitrary discrete Morse functions as out- 


4. 


5. 


6 . 


8 . 


lined in section 4 further strengthens our argument about simplicity of 


computing optimal discrete Morse functions. 


6 Concluding Remarks 

In this work, we provide a nearly linear time algorithm for computing ho¬ 
mology (with arbitrary coefficients) on 2-manifolds - the first such algorithm. 
This is particularly useful to compute homology of 2-manifolds that may 
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have torsion elements. The design involves the introduction and usage of 
an elementary simple homotopy construct that we call expansion frames. 
Having designed the optimal Morse function in linear time, we use a dynamic 
programming based pseudo-linear time boundary operator algorithm for com¬ 
puting the Morse boundary operator. Assuming the sum of Betti numbers 
is a small constant compared to the size of the complex, the Smith Normal 
Form is applied to a very smal input, giving us near-linearity. Finally, using 
the notion of expansion frames, we prove an unexpected result in discrete 
Morse theory: Start with an arbitrary DGVF on a 2-manifold and one may 
obtain an optimal DGVF merely by application of critical cell cancellations. 
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Appendix 

7 Elementary Algebraic Topology 

Definition 36 (Simplicial Complex). A simplicial complex 1C is a set of 
vertices and a collection C of subsets of vertices called faces. All faces 
satisfy the following property: The subset of a face is also a face. (i.e. 
B £ C,A £ B A £ C). Maximal faces w.r.t. inclusion are known as 
facets. The dimension of a face B is defined to be \B\ — 1. The dimension 
of the simplicial complex itself is the maximum over the dimension of its 
faces. 
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Figure 2: 2-manifolds 
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Definition 37 (Open Cell). An n-dimensional open cell is a topological space 
that is homeomorphic to an open ball. 

Definition 38 (Cell Complex). A hausdorff topological space X is called a 
finite cell complex if 

1. X is a disjoint union of open cells {Df} where Df is an open n-cell. 
(i ^ I where I is the indexing set.) 

2. For each open cell D'f there is a map —)■ X such that 

restricted to the interior of the closed ball defines a homeomorphism 
to Df and such that is contained in the {n — l)-skeleton of 

X. (The k-skeleton of X is the union of all open cells Di of dimension 
r < k). 

3. Finally, a set a is closed in X if and only if aCi Dj is closed in Dj 
for each cell . Note that Dj^ = (jfif{B'^). 

A cell complex is said to be regular if each is a homeomorphism and if it 
sends S'^~^ to a union of cells in the (n — l)-skeleton of X. 

In lay man terms, to construct a cell complex you start with points "D?, 
then glue on lines Vj to then glue discs Vf to Vj and T>^ and so on. 
Therefore a cell complex is a topological space constructed from a union of 
objects called cells, which are balls of some dimension, glued together on 
boundaries. Cell complexes are the most convenient object to do Algebraic 
Topology. But to simplify the discussion, we will instead provide a basic 
presentation of simplicial homology. 

Notation 39. Boundary & Coboundary of a simplex a: We define the 
boundary bd{a) and respectively coboundary Q a of a simplex as 
3 cj = {t\t < a} 5 a = {p\a < p} 

Homology groups are the most important and general topological in¬ 
variants of simplicial and cubical complexes, that are also computationally 
feasible. At the heart of it. Algebraic Topology is essentially the use of 
Linear Algebra to compute combinatorial topological invariants of a give 
space. Given a simplicial complex IT, can define simplicial g-chains, which 
are formal sums of g-simplices X]se5 where the Oj are integer coefficients. 
The abelian group of sums of fc-simplices under addition is called the Chain 
Group and denoted by Cq{W, Z). The n-simplex A = {uq, vi, - ■ ■ , Unjwith 
standard orientation is denoted -|- [uq, ui, • • • , Vn]- Consider the permutation 
group of n-letters on the vertices of A. The set of permutations fall into 2 
equivalence classes: even permutations and odd permutations. The set of 
even permutations induce the positive orientation -|- [uq, ui, • • • , Vn\ whereas 
the set of odd permutations induce the negative orientation — [uq, ui, • • • , Un]. 
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Figure 3: Dim 1 Boundary Operator 



Figure 4: Dim II Boundary Operator 

For each integer g, Cq{W) is the free abelian group generated by the set 
of oriented q-simplices of W. Let Wq be the total number of g—dimensional 
simplices for simplicial complex W. Then, one can show that Cq = 

The boundary map dq is defined to be the linear transformations dq : 

Cq ^ C,_l. 

Examples of such operations are given in Fig.ES and Fig.E4. 

This map gives rise to a chain complex: a sequence of vector spaces and 
linear transformations: 

0 ^ ^ Cn-l ... ^ Cq+l{W) ^ Cq{W) \ ... ^ Cl{W) % Cq{W) ^ 0. 

It can easily be proved that that for any integer g, 

dq o dq+i = 0. 

In general, a chain complex C* = {Cq,d} is precisely this : a sequence of 
abelian groups {Cq) connected by an operator dq : Cq ^ Cq-i that satisfies 
do d = 0. 
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V 2 

V 2 



• 

• 

^0 — 


Vo 

• 


Vl Vl 

Vl 

[ vo , Vi , V 2 ] 

.^[U0,ui] + [ vi , V2 ] 

Vl — Vo + V 2 — Vl 


+ b2,'f^0] 

+U0 - V 2 = 0 


Figure 5: dd = Ogives us a chain complex. 


If one defines 

Zq = kei dq and Bq = im9q+i, 

then it follows that Bq C Zq. Elements of Zq = kerclg are called cycles, and 
elements of Bq = im9g_|_i are called boundaries. Likewise, Zq = kerSg is 
called the q—th Cycle Group and Bq = imSg+i is called the q — th . Boundary 
Group. Then the homology group Hq measures the equivalence class of cycles 
by quotient-ing out the boundaries i.e. this construction measures how far 
the sequence is from being exact. 

The g-dimensional homology of W , denoted Hq{W) is the quotient vector 
space. 


Hq{W) 


Zq{W) 

Bq{W) 


and the ( 7 -th Betti number of W is its dimension: 


/3q = dim Hq = dim Zq — dim Bq 


8 Morse Homology 

Let F be a Discrete Morse function defined on simplicial complex W. Let 
Cq{W, Z) denote the space of ( 7 -simplicial chains, and Mq which is a subset 
of Cq{W, Z) denote the span of the critical g-simplices. Let denote the 
space of Morse chains. Let Cq denote the number of critical g-simplices. Then 
we have, M.q = Z'^?. 
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Theorem 40 (Forman | 6 ]). There exist boundary maps dq : Mq —?■ Mq-i, 
for each q, which satisfy 

dq o dq+i = 0 . 

and such that the resulting differential complex 


t 4 dn 1 j dn—l 

Mn -^Mn-1 




82 , j d\ ,4 

Ml —^ Ml) 


0 -^ jvin -^ J^ln-1 -^ • • • -^ 

calculates the homology ofW. i.e. if we go with the natural definition, 

imOq-il 


Then for each q, we have Hq{M,d) = HqiyV^Tj). 

Theorem 41 (Boundary Operator Computation. Forman |6]). Consider an 
oriented simplicial complex. Then for any critical (p+1 )-simplex ft set: 


dp = E Pah « 

critical ck(p) 

Pah = E ^^(7) 

■y&r{h,a) 

where r(P,a) is the set of discrete gradient paths which go from a face 
in d P to a. The multiplicity Npy) of any gradient path 7 is equal to ±1 
depending on whether given 7 the orientation on P induces the chosen orien¬ 
tation on a or the opposite orientation. With the boundary operator above, 
the complex computes the homology of complex K. 


Theorem 42 (Forman | 6 ]). If a < b, are real numbers, such that fa,b] 
contains no critical values of Morse function {, then the sublevel set M(b) is 
homotopy equivalent to the sublevel set M{a). 


Theorem 43 (Forman|6j). Suppose is a critical cell of index p with 
/(ct) G [a, 5] and f~^{a,b) contains no other critical points. Then M{b) is 
homotopy equivalent to 

M{a) IJ eP 

< 

where e^ denotes a p-dimensional cell with boundary e^. 


In Thm,43 Forman’s establishes the existence of a cell complex (let 
us call it the Morse Smale Complex) that is homotopy equivalent to the 
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Figure 6: Frame Expansion: Example 1. Part L 


original complex. For proof details please refer to FormanjH]. The boundary 
operator in Thm,41 for the chain complex construction (referred to as the 
Morse complex) tells us how to use the new CW complex that is built in 
construction described in proof of Thmj4^ Note that the Morse complex 
itself is a chain complex and not a CW complex. But, the chain complex 
construction (referred to as the Morse complex) tells us that both these 
constructions have identical homology. 


9 Extra Figures 

10 Detailed Pseudocode 
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Figure 9: Frame Expansion: Example 2: Part II 



Figure 10: Frame Expansion: Example 2: Part III 
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Figure 11: Frame Expansion: Example 3: Part 1 
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Figure 12: Frame Expansion: Exampe 3: Part 2 
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Algorithm 4 Homology() 

1: procedure CALcHoMOLOGY(At, A;) 

2: y := mainFrame(At); 

3: Ac := calcbdryOp(Al,^, y); 

4: ;= SmithNormalForm(Ac, A); 

5: end procedure 


6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


procedure CALcBDRYOp(At, 1^) 
topologicalSort(^, y, £, 'ASCENDING'); 
for 1 < i < |£|; a := C[i] do 

if dim O' = 0 & cj-pair = NIL then 
Act = 0; 

else if fj ^ /3 & cj-pair = /3 then 

Act =< 9/3, fj > xA/3; 

else 

Let Tj -< cr be the set of regular cells incident on a s.t. {Ti,a) ^ Y] 
Let Q!i -< cr be the set of critical cells incident on rr; 

Act = ^ < 9a, r* > + X) cri x < da, a* >; 

end if 

if cr-pair = NIL 8^ rr-revPair = NIL then 
AcCJ := Act; 

end if 
end for 

end procedure 


39 





Algorithm 5 Algorithm FrameFlow() 

1: procedure findCoBdry(A1'^,;B‘^) 

2; Go through list If has only one face, add to B 

3; return B\ 

4: end procedure 

5: procedure ADDPAiRToVECTORFiELD(r, -d, 'f ,M,B‘^, d) 

6: if r-revPair = NIL then 
7: delete(r, delete(i9, 

8: if r G B then delete(r, ,6); end if 

9: nQ(Q, r); •d- pair := r; r- revPair := -i?; 

10; y --y + 

11: end if 

12: end procedure 

13: procedure frameFlow(A 1, d, A) 

14: if (tz7 = dQ(,B'^)) = NIL then 
15; I VD ;= dQ(At'^); 

16: end if 
17: repeat 
18: fF := 'T'; 

19: repeat 

20: A := faces(cL7) — tu-revPair; v := bdry(A); 

21: iffF='T' NIL then 

22 ; I addPairToVectorField(ti7, R, A, At, d); 

23: else 

24: I delete(cr 7 , At*^); 

25: end if 

26: fF:='F'; 

27: for 1 < i < |A|; do 

28; d := A[i]; fx := cofaces(i?) — w\ 

29: addPairToVectorField(ii, i?, A, At, d); 

30: end for 

31: until {w = dQ(Q)) / NIL 

32: until (tj7 = dQ(,B'^)) ^ NIL 

33: end procedure 

34; procedure MAiNFRAME(At, A) 

35; B"^ := findCoBdry (At^); 

36; frameFlow (At,2,,6^,A); 

37: {At^[l : numEars], numEars} := earDecompose(At); 

38: for 1 < i < numEars do 40 

39: := findCoBdry(At^[i]); 

40: frameFlow (At^[i],l,;B^[i],A); 

41: end for 
42: return A.; 

43: end procedure 







